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Abstract 

The trace formula for the density of single-particle levels in the two-dimensional radial power- 
law potentials, which nicely approximate the radial dependence of the Woods-Saxon potential and 
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■rj" , quantum spectra in a bound region, was derived by the improved stationary phase method. The 



.^ ' specific analytical results are obtained for the powers 4 and 6. The enhancement phenomena near 

o ■ 

the bifurcations of periodic orbits are found to be significant for the description of the fine shell 



m 



structure. It is shown that the semiclassical trace formulas for the shell corrections to the level 



S^ ] density and energy reproduce the quantum results with good accuracy through all the bifurca- 

H ■ 
_Cp_. tion (symmetry breaking) catastrophe points, where the standard stationary-phase method breaks 

down. Various limits (including the harmonic oscillator and the spherical billiard) are obtained 

from the same analytical trace formula. 
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I. INTRODUCTION 



According to the shell-correction method (SCM) [l|, |2|, the oscillating part of the total 
energy of a finite fermion system, the so-called shell-correction energy 5U , is associated with 
an inhomogeneity of the single-particle energy levels near the Fermi surface. Depending on 
the level density at the Fermi energy - and with it the shell-correction energy 5U - be- 
ing a maximum or a minimum, the nucleus is particularly unstable or stable, respectively. 
This situation varies with particle numbers and nuclear-potential parameters like deforma- 
tion ones. In consequence, the stability of nuclei depend strongly on particle numbers and 
parameters of the potential. 

A semiclassical theory of shell effects, the so-called periodic orbit theory (POT) 3|-[7|], 
was used for a deeper understanding, based on classical pictures, of the origin of nuclear 
shell structure and its relation to a possible chaotic nature of the dynamics of nucleons. 
This theory provides us with a nice tool for answering, sometimes even analytically, the 
fundamental questions concerning the physical origins of the double-humped fission barrier 
and, in particular, of the existence of the isomer minimum 8l-ll3l| (see Refs. 12|, [iSj for the 
introductory review and text book). Some applications of the POT to nuclear deformation 
energies were presented and discussed for the cavity (billiard-like) potentials with sharp 
edges in relation to the bifurcations of periodic orbits with the pronounced shell effects. 

In the way to more realistic semiclassical calculations, it is important to account for a 
diffuseness of the nuclear edge. It is known that the central part of the realistic effective 
mean-filed potential for nuclei or metallic clusters are described by the Woods-Saxon (WS) 
potential Vws('^) \lm- The idea of Refs. [l5|, |l6| is that the WS potential is nicely approxi- 
mated (up to a constant shift) by much a simpler power-law potential which is proportional 
to a power of the radial coordinate r°. The approximate equality 



Vws(r)^VWs(0) + W^or" 



:i.ii 



holds up to around the Fermi energy with a suitable choice of the parameters Wq and a. 
In the case of the spatial dimension T) = 2, one can use Eq. 01.11) for a realistic potential 



of electrons in a circular quantum dot [12|, Il3|, [IT] . We shall derive first the generic trace 
formula for this radial power-law (RPL) potential in the case of two dimensions, and then 
discuss its well known limits to the harmonic oscillator and cavity (billiard) potentials [l3|. 



The manuscript is organized as follows. In Sec. |TT]the classical dynamics is specified for the 
RPL potentials. The outlook of the general phase-space trace formula for arbitrary potentials 
is presented in Sec. Illli The derivations of the trace formulas for the RPL potentials in two 
dimensions is given in Sec. IIVI Section |V] is devoted to the comparison of the semiclassical 
calculations for the oscillating level density, and the shell-correction energy with quantum 
results. The paper is summarized in Sec. I VII Some details of our POT calculations, in 
particular full analytical derivations at the powers a = 4 (see also Ref. [l8|) and 6 for all 
periodic orbits (POs) and those at arbitrary a for the diameter and circle orbits, are given 
in the Appendixes RHEl 

II. CLASSICAL DYNAMICS AND BIFURCATIONS 

The radial power-law (RPL) potential model is described by the Hamiltonian 

where m is the mass of the particle, -Rq ^"^^ Eq are introduced as constants, having the 
dimension of length and energy, respectively, and are related with Wq in Eq. (II. ip by Wq = 
Eq/Rq. (Thus, Eq can be arbitrary, and it is convenient to take Eq = h'^/niRl in quantum 
mechanical formulation.) This Hamiltonian includes the limits of the harmonic oscillator 
[a = 2) and the cavity {a — )■ oo); realistic nuclear potentials with steep but smooth surfaces 
correspond to values in the range 2 < a < oo. The advantage of this potential is that it 
is a homogeneous function of the coordinates, so that the classical equations of motion are 
invariant under the scale transformations: 

with E ^ sE . (2.2) 

Therefore, one only has to solve the classical dynamics once at a fixed energy, e.g., E = Eq 
(s = 1); the results for all other energies E are then simply given by the scale transformations 
with s = E/Eq hy definition in the last equation of Eq. (12. 2p . This highly simplifies 



the POT analysis 15|, ll9|. Note that the definition (12. ip can also be generalized to include 



deformations (see, e.g., Ref. |l6l. Il9j) 



As the RPL Hamiltonian H is spherically symmetric, it can be written explicitly in the 
two-dimensional (2D) spherical canonical phase-space variables {r, y9;pr,Pv'}' 

where ip is the azimuthal angle (a cyclic variable), p;^ = L is the angular momentum, and 
the radial momentum p^. is given by 



Pr{r,L)- ^p^{r) ^2 
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Pir) = \hm 


E- 


-Eo 





(2.4) 

The classical trajectory (CT) r(t) can be easily found by integrating the radial equation of 
motion r = Pr/m with Eq. 02 ■4p . Transforming the spherical canonical variables into the 
action- angle ones, for the actions /,., I^p one has 

Prdr =Ir{E,L), (2.5) 



vr ,„ 

" ' n 



1 r^TT 



I^ = — p^dip = L , (2.6) 

^TT Jq 

where r^^^ and Tmax are the turning points which are the two real (positive) solutions of the 
equation p^(r, L) = 0. 

The definition (12. ip can be used in arbitrary spatial dimensions, as long as r is the 
corresponding radial variable. In practice, we are interested only in the 2D and 3D cases. 
The spherical 3D and the circular 2D potential models have common sets of periodic orbits, 
see Fig. [H For a > 2, POs with the highest degeneracy (/C = 1 and 3 in the 2D and 3D 
cases, respectively) are specified by three integers and labeled as M{nr,n^), where Ur and 
n^ are mutually commensurable integers. They represent the number of oscillations in the 
radial direction (n^), and the number of rotations around the origin (n^) for one period, 
and M is the repetition number. For the isotropic harmonic oscillator [a = 2), all the 
classical orbits are periodic ones with (degenerate) ellipse shapes. By slightly varying a 
away from 2, the specific diameter and circle orbits appear separately, and they remain as 
the shortest POs with the corresponding degeneracies /C = 1 and 0. With increasing a, the 
circle orbit and its repetitions cause successive bifurcations generating various new periodic 
orbits {nr,n^}, rir > 2n<p. Fig. [1] shows some of the shortest POs M{nr,n^)- The shortest 
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FIG. 1. The scaled periods t^q of some short periodic orbits, plotted as functions of the power 
parameter a in dimensionless units m = K^^ = Eq = \ (Appendix B). Thin solid curves are the 
circle orbits MC ^ dashed curves are the diameters Af (2, 1), and thick solid curves are the polygon- 
like orbits M{nr,n^) {rir > 2n^), which bifurcate from the circle orbits MC at the bifurcation 
points indicated by open circles. 

PO is the diameter which has the degeneracy /C = 1 in the 2D problem at a > 2. Other 
polygon-like orbits have /C = 1 at a > abif, where abif is a bifurcation value (see its specific 
expression below). The circle orbit having maximum angular momentum is isolated (/C = 0) 
for the 2D system (except for the bifurcation points). 

For the frequencies of the radial and angular motion of particle, one finds 



Ur 



dir \dE' ' - 



-1 



OH f dlr\ ^ ( vn 



L 

dH 



dr 



w^ 



{dlr/dL)^ 



dL {dlr/dE)^ ' ^^■'^^ 

where Pr{r,L) is the radial momentum, given in Eq. (12. 4p . and J^ = Ij.{E,L) [Eq. (12. 5p ] is 
identical to the energy surface H{Ir, L) = E. Thus, the PO condition is written as 



L 

71 



dr 



n„ 



2 r n • (2.8) 

The energy surface /^ = Ir(E,L) is simplified to a function of the only one variable L 
[see Eq. (12. 5p ]. The solutions to the PO equation (12. Sp . L* = L*{nr,n^), for the given co- 
primitive integers n^ and rir define the one-parametric families /C = 1 of orbits M{nr,n^) 



because L is the single- valued integral of motion, which is only 1, besides of the energy E. 
[Notice that another single- valued integral of motion, Q = sin{nr0ip — n^Or) , appears under 
the PO condition fl2.8p {Or and 6*^ are the angles of the action-angle variables). However, 
for these RPL potentials, no additional degeneracy arises since the rationality of frequencies 
occurs only at the discrete values of angular momenta L*; see Refs. p, |8|.] The azimuthal 
angle (f can be taken, for instance, as a parameter of the orbit of such a family. 

According to the limit f{L) — )■ 1/2 at L — )■ 0, one has the diameter orbits M(2, 1) as the 
specific one-parametric (/C = 1) families related to the solution L = of Eq. (12. 8p . The other 
specific solutions are the isolated (/C = 0) circle orbits MC by which we represent the M-th 
repetition of the primitive circle orbit C. The radius r^ of the circle orbit is determined 
by the system of equations, namely r^^^ = rmax = r^ for equivalence of the turning points 
being, in turn, the solutions of the equation p1{E,L) = (see Appendix R|) . Thus, one has 
also the angular momentum of the circle orbit Lc = r^p{r^). As seen obviously from the 
condition of the real radial momentum pr [Eq. ( 12. 4^ ]. this Lc is the maximal value of the 
angular momentum L, i.e., < |L| < Lc ■ 

As shown in Appendix lAj for the calculation of the Gutzwiller stability factor Fmc of 
the circle orbit MC in the radial direction, one obtains 

'ttMHc' 



Fmc = 2 - Tt{Mc)^'^ = 4sin2 
= 4 sin^ TrMV2 + a 



(2.9) 



where Qcl'^c ^^ '^^ ratio of the radial-to-angular frequencies. 



V2T^, o;^ = a;^(L = Lc) = — V • (2-10) 






The factor Fmc becomes zero at the bifurcation points abif by the definition of the stability 
matrix Tr(A^c')^^ = 2; see Eq. (12.91) . According to the first equation in Eq. (12.101) . one finds 



nt 



«bif = ^-2 (2.11) 

for the POs, Qc/^c ~ f^r/n^. The PO family M{nr,n^), which corresponds to the so- 
lutions L* < Lc of the PO equation (12. 8p . exists for all a > aui- There is one specific 
bifurcation point a = 2 in the spherical harmonic oscillator (HO) limit with the frequency 
u^ = ^J2EQ/{mB^)^ where one has the two-parametric families at any L within a contin- 
uum < L < E /uju,- In this limit the above specified circle and diameter orbits belong 



to a two-parametric family mentioned above. For the circular billiard limit a — )■ oo, the 

isolated circle orbit (/C = 0) disappears degenerating into the billiard boundary r^ ^ Rq, 

Lc — !■ v2rnE Rq [see the limit a — )■ oo in Eq. (1A2I) for r^]. 

Another key quantity in the POT is the curvature K of the energy surface Ir = IriE, L) 

given by 

d^Ir{E,L) df{L) 

where f{L) is the ratio of frequencies defined in Eq. fl2.8p . As shown below, the curvature 
f l2.12p and Gutzwiller factor (12. 9p . related both to the stability of the orbit, are the key 
quantities for calculations of the magnitude of the PO contributions into the semiclassical 
level density. 

III. PHASE-SPACE TRACE FORMULA 



The level density g{E) for the Hami 
space trace formula (in T) dimensions) 



tonian ij"(r, p) can be obtained by using the phase- 



lol, 
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20 



23|: 



9UE) = j^ Re ^ I dr' J dp" 6 {E - H{v'\ p")) 

X \Jct{v'l.v1)\"^ exp r^^cT - ^|/icx) • (3-1) 

The sum is taken over all discrete CT manifolds for a particle moving between the initial 
r', p'; and the final r", p" points with a given energy E. Any CT can be uniquely specified 
by fixing, for instance, the initial condition r', and the final momentum p" for a given time 
tp^ of the motion along the CT. For the action phase ^ct in exponent of (13. ip . one has 



^CT = 5cT(p',p",tcT) + (p"-p')-r' 

= 5cT(r',r",E)-p"-(r"-r'), (3.2) 



where iSctIp', p"5^ct) ^'^'^ •S'ctIi"', r", -E) are the actions in the momentum and coordinate 
representations, respectively, 

r" 

Sct{v\ p", ^ct) = - / dp ■ r(p) , 

r" 

ScT{r',T",E) = j ^ dr-p(r). (3.3) 



In Eq. (13. ip . J7ct(Pj_, P±) is the Jacobian for the transformation of the initial perpendicular- 



■ov phase related to the number of 

2a. 



to-CT momentum p^ to the final one p" . fi^rj, is the Mas^ 
conjugate (turning and caustics) points along the CT [24 

For calculations of the trace integral (13. ip by the stationary phase method (SPM), one 
may write the stationary phase conditions in both p" and r' variables. According to the 
definitions (13. 2p and (13. 3p . the stationary phase condition for the p" variable is that for the 
CT to be closed in the coordinate space: 



dp" J -- '") =°- ^^-^^ 

Here and below, the star means that any quantities (as in the parentheses) are taken at the 
stationary point; in particular, p" = p"* in the variable p". In the next integration over r' 
by the SPM, we use the Legendre transformation (13. 2p . and the closing condition (13. 4p which 
lead to <?CT = ScT{'r',r",E). Therefore, the stationary phase condition for this integration 
over the spatial coordinates r = r' = r" can be written in the following forms: 



dr' dr" J \ dr' dr" 

= - (p' - P")* = , (3.5) 

where the star means r' = r" = r* along with p' = p" = p*. Therefore, one has the 
condition for a CT to be closed also in the momentum space, and then, the stationary phase 
conditions are equivalent to the periodic- orbit equations (13. 4p and (13. 5p . After appl yin g these 



two conditions we may arrive at the trace formula in terms of the sum over POs 12|, |l3| 



IV. TRACE FORMULA FOR THE 2D POWER-LAW POTENTIALS 

As showed in Sec. Ull there are two kind of POs in the 2D circular potential; i.e., one- 
parametric families (/C = 1, including the diameters), and isolated (/C = 0) circle POs; for 
which we shall derive the trace formulas separately in the subsections A and B, respectively. 
In these derivations for both kinds of POs, it is useful to take the canonical integration 
variables with the angular momentum L as one of the integration variables in the phase-space 
trace formula (13. ip . As shown in the subsection A for the contributions of the one-parametric 
(/C = 1) families of the diameters and polygon-like POs, the integration over the cyclic 
angle variable ip, conjugate to L, is easily carried out (it simply gives 27r). The remaining 



integration over L is done within the simplest improved SPM (ISPM [ll|, |12|, |22|, |23|) with 
the finite hmits —Lc < L < Lc- Thus, we shall obtain the ISPM trace formula in terms of 
error functions 12|. Extending the integration limits to ±00 far from the critical symmetry- 
breaking points, we will arrive at the asymptotic Berry- Tabor trace formula |6| for the 
contributions of P (polygon-like) and D (diameter) POs. In the next subsection B, the 
contribution of the isolated (/C = 0) circle orbits is derived by adopting the ISPM for 
two integrations over suitable (radial phase-space) variables. It yields a product of two error 
functions. Extension of the finite integration limits in both error functions to infinity (i.e., to 
±00 after their transformation to the corresponding Fresnel integrals of positive or negative 
real arguments) leads to__the standard SPM (SSPM) Gutzwiller trace formula |3[ for the 
considered isolated POs 

A. One-parametric orbit families (/C = 1) 

For the contribution of the one-parametric families of the maximal degeneracy /C = 1 into 
the phase-space trace formula (13. ip , it is useful to transform canonically the usual Cartesian 
phase-space variables {p; r} to the action-angle ones {I; ©}, specified in the spherical action- 
angle variables as = {OriOip = v'}; I = {Irilip = L}. The Hamiltonian if, action phase 
<PcT, and other related quantities of the integrand in Eq. (13. ip are functions only of action 
variables [H = H(I) = H{Ir, I^) = H{Ir, L)], i.e., they are independent of the angle variables 
0. Therefore, performing easily the integrations over these angle variables 0, one has the 
factor 27r per each integration. Taking then the integral over /^ exactly by using the energy 
conserving (5-f unction, one obtains 






X exp l—M K h{E, L) + n^L] - — /i,,„^ „^ j . (4.1) 

Here, the phase (13. 2p is expressed in terms of the corresponding action-angle variables 
through the actions (13.30 in the considered mixed representation, 

^CT = 27rMKJ,(E,L)+n^L], (4.2) 

Ur and n^ are positive co-primitive integers, M is a nonzero integer, uj^ is the radial frequency 
in Eq. (12. 7p . We took also into account explicitly that the actions I (or {/^, L}) are constants 

9 



of motion for the spherical integrable Hamiltonian, omitting in I the upper indexes, related 
to their initial (prime) and final (double prime) values of Eq. (13. ip . The integration limits 
in Eq. (14. ip for L are —Lq < L < Lq-, where Lc is the maximum value corresponding to 
the circle orbit. All quantities in the integrand are taken at the energy surface /^ = It{E, L) 
[Eq. (12.51) ]. Thus, Eq. 04.11) is similar to the semiclassical Poisson summation trace formula 
which can be obtained directly by using the EBK quantization rules [6|, 113|] for the case 
of the spherical symmetry of the Hamiltonian. Note that formally, before taking the trace 
integral over the angular momentum L by the SPM in Eq. (14.10 . one can consider positive 
and negative M, as those related to the two opposite directions of motion along a CT. 
They give, of course, equivalent contributions into the trace formula, due to a time-reversal 
symmetry of the Hamiltonian [or one may write simply the additional factor 2 in Eq. (14. 1|) . 
but with a further summation over only positive integers M]. However, in contrast to the 



standard Poisson summation trace formula 



13|, there is no zero values of the integers in 



Eq. ( 14. ip . n^/rir > 0. In the derivations of Eq. ( 14. ip from Eq. ( 13. ip . we essentially used 
that for families of the maximal degeneracy /C = 1, the generating function ^ct [Eq. ( 14. 2p ] 
becomes independent of the angle variables for the integrable Hamiltonian. Notice that in 
these derivations, the SPM conditions (13. 4p and ( 13. 5 p were satisfied identically within the 
continuum of the stationary points < ^p,Or < 27r, which form CTs, but they are not yet 
POs generally speaking for arbitrary angular momentum L. (Exceptions are the cases of the 
complete degeneracy like the spherical HO; see below.) The integration range in Eq. (14. ip 
taken from the minimum, L_ = 0, to the maximum, L+, value (for anticlockwise motion, 
for instance) covers the contributions of a whole manifold of closed and unclosed CTs of the 
tori in the phase space at the energy surface around the stationary point, L = L*, which 



corresponds to the PO 23|. We shall specify the integration limits L+ for the contribution 
of the (/C = 1) diameter families M{nj. = 2,n^ = 1) into Eq. (14. ip in Appendix [Dl 

We apply then the stationary phase condition with respect to the variable L for the 
exponent phase <Pct [Eq- (14. 2p ] in the integrand of Eq. ( 14. ip . 

{d^c^/dLf = 0, (4.3) 



which is exactly the resonance condition ( 12. 8p . This condition determines the stationary 
phase point, L = L* = Lpo, related to the POs M{nr.,n^) of the families. Expanding now 
the exponent phase <Pct [Eq- (14. 2p ] in the variable L up to the second order, and assuming 

10 



that there is no singularities in the curvature f l2.12p for the contribution of all /C = 1 families, 
one has 

^CT = SpoiE) + ^ jS(^ - L*f + ■ ■ • , (4.4) 

where Spo{E) is the action along one of the isolated PO families, determined by Eq. f l2.8p . 

^Po(^) = 27rM [ur Ir{E, L*) + n^ L*] . (4.5) 

In this equation, M is the number of repetitions along the PO = M{nr, n^), Ir{E, L) is the 
energy surface [Eq. (12. 5p ]. L = L*{nr,n^) is the solution of the PO equations (12. 8p or (14. 3p . 
The Jacobian JpQ in Eq. (14. 4 p measures the stability of the PO with respect to the variation 
of the angular momentum L at the energy surface, 

4o = (^1^) = 2vrMn,i^po, (4.6) 



where Kpo is the curvature (I2.12p . (IBSp of the energy surface J^ = Ir{E, L) at L = L* = LpQ. 
For the sake of simplicity, we shall discuss the simplest leading ISPM with the second 
order expansion (14. 4 p of the action phase, and with the zero order one of the pre-exponent 
factor over the L variable in L — L* in Eq. (14. ip . Substituting now the expansion (14. 4p 
into Eq. (14. ip . one takes there the pre-exponential factor off the integral at L = L*. Thus, 
applying Eq. (14. 4p . we are left with the integral over L of a Gaussian type integrand within 
the finite limits mentioned above for contributions of the one-parametric polygon-like and 
diameter families, including the contribution of boundaries for < n^/nj, < 1/2. Taking 
this integral over L within the finite limits, one obtains the ISPM trace formula, Sg^'^\E), 
for contributions of the one-parametric (/C = 1) orbits. 



^'\E) = ReJ24'kE) 



^(1 ,^j - ... 

PO 



X exp 



z 7r 

- Spo{E) - i-apo - i(f)d 



(4.8) 



The sum is taken over the discrete families of the PO M{nr,n^) with rir > 2n^, M > 1 in 
the 2D RPL potential. Spo(E) is the action fl4.5p a long these POs. For the amplitudes ApQ 
as for the /C = 1 families in the elliptic billiard |22|, and the integrable Henon-Heiles (IHH) 



potentials |23|, one finds 



4o = , ; "'^ erf (Zpo, Z+o) (4.9) 
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with the period 

Tpo = ^^ = ^^ . (4.10) 

In Eq. f l4.9p . K-pQ, is the curvature of the energy surface /^ = Ir{E,L) [see Eqs. f l2.12p . flBSp ], 
The generahzed complex error function in Eq. (14. 9 p is introduced by 

2 r _ 2 

erf (u, t;) = —^ / dz e "" = erf(w) — erf (u) (4-11) 



with the standard error functions, erf (2;), of the complex arguments z. These arguments are 
specified by 



Z^Q = ^/-i7iMnr\Kpo\/h (L± - Lpo) , (4.12) 

L„ = and L+ = Lc for all /C = 1 polygon-like PO families (besides of the diameters, 
see below). For simplicity, the finite integration interval of the angular momenta was split 
into two parts, —Lc < -^ < and < L < Lc, where Lc is the angular momentum 
of a circle orbit, as mentioned above. There are the symmetric stationary points, ±|i^*|, 
related to the anticlockwise and clockwise motions of the particle along the PO in two 
these phase-space parts. They give equivalent contributions to the amplitude, due to the 
independence of the Hamiltonian of time. Thus, we have reduced the integration region 
to < L < Lc, accounting for this time- reversibility symmetry simply by the factor 2 in 
Eq. (14. 9 p (exclusion is the diameters, for which one has one stationary point L* = 0, and 
therefore, the time-reversibility degeneracy, taken into account automatically by the limits 
of the error functions, is one). For all the polygon-like and diameter POs (n^ > 2), we also 
found L_ = for the minimum value of the angular momentum L. 

For the Maslov index of the considered K, = 1 PO families and the constant phase 0^ in 
Eq. (14. Sp . one obtains 

a^il = 2Mnr, ^d = -7r/4 . (4.13) 

The Maslov index OpQ is determined in terms of the number of turning and caustic points by 



the Maslov&Fedoryuk catastrophe theory, see Refs. |4j. |23|-|26I| . Note that for the potentials 



with smooth edges, one has another expression for the Maslov index ap„ than for the circular 



billiard 



13 



271]. Note also that the total Maslov phase, defined as a sum of the asymptotic 

, depends on the energy 
22|]. This total Maslov 



part (I4.13P and the argument of the complex density amplitude f | 4.9l 



E and parameter a of the RPL potential [Eq. (11.10 : see Refs. H, 

phase is changed through the bifurcation points smoothly, due to the phase of the complex 

error function in the amplitude (14. 9 p in Eq. (14.80 . 

12 



For the stationary point L* far from the ends of the physical integration interval, one can 
extend the integration range to the infinity from — oo to oo (in the case of diameters from 
zero to oo). We then arrive asymptotically at the Berry&Tabor result |6| for the contribution 
of all /C = 1 families f l4.8p with the following amplitude: 



^PO 



Axm ^ 



-ni 



where d-po takes into account the discrete degeneracy, dp q = 1 for diameters M(2, 1) [n. 



(4.14) 



2n^), and 2 for all other (polygon-like) POs (n^ > 2n^) 13|. In the circular billiard limit 

(a — )■ oo), the action Spo{E) is given by 



Spo{E) = pCpo, p = v2mE 



(4.15) 



with Cpo being the length of the PO. For the curvature Kpo [see Eqs. (I2.12p and (IBSp ]. one 
can asymptotically {a — > oo) obtain 

^po = „ . ] ^^ • (4.16) 

Substituting the quantities (I4.15p . (I4.16P and (I4.13P into Eq s. (14. 8p and (I4.14p . one obtains 



the well known trace formula for the circular billiard [13|, l27|. Note that the amplitude 
(14. 9 p of the solution (14. 8 p is regular at the bifurcations which are the boundary points 
L = L* = Lq of the action (L) part of the tori as in the elliptic billiard 22| . 

Our SSPM result (I4.14p coincides with the Berry and Tabor trace formula J6[ , as adopted 
to the 2D spherically-symmetric Hamiltonians by using the simplest expansions of the action 
phase and amplitude near the stationary point (see above), instead of a more general but 



more complicated mapping procedure; see more comments in Ref. 22|. The essential differ- 



n 



ence from the Berry&Tabor theory [6| is that Eq. (14. 8 p is one of the terms of solutions of the 
breaking-of-symmetry problem for the highest degenerate orbits, such as the one-parametric 



families in the IHH potential, or the elliptic and hyperbolic orbits in the elliptic bil 



(see also Refs. 



12 



23|). Within the SPM of the extended Gutzwiller approach 



iardj22| 

yy, 



we have to derive separately the contributions of the other orbits as the circle /C = POs 
in the RPL potentials beyond the semiclassical Poisson summation-like trace formula (14. ip 
(with the restrictions to the range of the Ur and n^ integer variables). Note that the ISPM 
trace formula (14. 8 p for the one-parametric families contains the end contributions related to 
the finite limits of integrations in the error functions. This trace formula is the contribution 
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of only such families, as essentially was used in its derivation from the trace formula (13. ip 
by taking 27r per the integral over the angle variable. Therefore, there is no contributions of 
the circle orbits in Eqs. (14.1 p and fl4.8p . These orbits correspond to the isolated stationary 
phase point, but the derivation of this section HV Al cannot be applied for their contributions 
within the ISPM. This is similar to that for the IHH potential 



Notice also that all roots L = L* = Lpo of the stationary phase equation (12 .Sp for /C = 1 
families PO = M{nr,n^) are in between the minimum value L = L* = for diameters, 
and a maximum one L = Lc, < Lpo < Lc (anticlockwise motion, for example). There 
is the additional isolated stationary-boundary point L* = Lc, which gives the separate 
contribution of the circle orbits into the phase space trace formula, see below. 

B. Circle orbits (/C = 0) 

In contrast to the derivations of contributions of the most degenerated (/C = 1) orbits, we 
take now into account existence of the isolated stationary point of the phase ^ct (13. 2p in the 
radial spherical phase-space variables r' * = r" * = r^, p'^* = p'^* = [Eqs. (13. 4p and (13. 5p ]. 
After the transformation of the integration variables in Eq. (13. ip to the spherical phase 
space coordinates {r',if';p'^,L}, it is convenient first to perform the exact integrations over 
L by using the energy conserving (5-function, and over the cycle azimuthal angle ip' leading 
simply to 27r as above {j dip' /uo^p = T<^,ct is the primitive rotation period). Thus, as in the 
subsection A, one finds 



9sci{E) = j^ ReJ2 jdr' jdp':.T^,cT 



X \JcT{p'r,Pr)\ exp 



1/2 



I IT 



(4.17) 



The stationary phase condition for the SPM integration over the radial momentum p'^ in 
Eq. fHTTl) writes 

= (r' - r") = . (4.18) 



dp'l 

The solution of this equation is the isolated stationary point p'^ = p'l* = P* = 0. The phase 
<PcT [Eq. (13. 2p ] is expanded in the momentum p'l near this point p^' * = in power series, 

^CT = ^CT + \jS^ {P'r - Plf + ■ ■ ■ , (4.19) 
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where the Jacobian J^rj. is given by 






2nMnrK ' * 
{dp'l/dLf 



The upper index star means again that the corresponding quantity is taken at the stationary 

point, p'^ = p'r* = 0. By using the 2nd order expansion of the exponent phase fl4.19p and 

taking the pre-exponent amphtude factor off the integral at this stationary point, one gets the 

internal integral over p" in terms of the error function as in the previous section. According 

to Eq. (13. 2p . with the radial-coordinate closing condition (14.18^ for the CTs, the short phase 

<Pct i^^ Eq. f l4.19p can be written in terms of the corresponding variables as ^^t — I > Pr dr. 

Taking then into account the CT closing condition f l4.18p . r' = r" = r, for the stationary 

phase equation in the integration over the radial r coordinate perpendicular to the circle 

orbit, one results in 

»-,»ax^•_(p,,_p,).^0, (4.21) 



dr" dr' 

Therefore, together with Eq. (14.18^ . one has the PO conditions related to the circular orbit 

r = r* = r^ and L = L* = Lq (see Appendix |X]). As usually within the SPM, we expand 

now the phase ^^ in the radial coordinate r near this r* = r^ 



'c 



^hT = MSc + ^ j£J, (r - rcf + ■ ■ ■ , (4.22) 



where Sc is the action along the primitive circle PO (C), 

^(r) _ f d ScT _^ ^ d 5'cT ^ d ScT 



Mc \ Q^/2 ^ dr'dr" Q^n'i 



MC 



dp'^ ^M + M^ . (4.23) 



dr' dr" dr" , ^^^ 



Using now the action phase expansion (I4.22p at the second order as the simplest ISPM 
approximation, and taking the pre-exponent amplitude factor at the isolated stationary 
point r = r(7 off the integral, one finally obtains 



^°^^^(E) = Re5^A£UE) 



M=\ 



X exp 



(4.24) 



The sum runs all circle orbits MC^ M = 1,2, ■■ ■ are positive integers. The action Sc{E) 
along the primitive C orbit is given by 

Sc{E)= (I p^d^ = 27iLc (4.25) 
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with Lc^ given explicitly in Eq. flA2p . In Eq. f l4.24p . (TjJ(j is the Maslov index determined by 
the number of caustic and turning points along the circle orbit, according to the Fedoryuk& 
Maslov catastrophe theory (see Refs. |23|-|26[|). 



"MC = 4Af, 



,(0, 



. 



(4.26) 



For the amplitudes A)^lfj[E) in Eq. 04.24p . one finds 



A 



(0) 
MC 



T, 



c 



eii{Z. 



p,MC^^p,MCJ ^''^^K^r,MC^ ■^r,Mc) ' 



(-) ^(+) 



(4.27) 



(4.28) 



where Tc is the period of the primitive (M = 1) orbit C 

Tc = 2iTmrl/Lc , 



see Eq. (]A2p for r^ and Lc- In Eq. (I4.27p . Fmc is the Gutzwiller stability factor [3[ of 
the circle orbits [Eq. (12. 9p ]. The arguments of the error functions in Eq. (I4.27P can be 
transformed to the following invariant form (see Appendix lEl): 



7(±) 

^p,MC 



h 



IT 



MVa + 2 



'C, 



Kc 



0, 



^c) 



Z, 



(±) 

MC 



iF. 



MC 



\ 47rMn(a + 2)3/2 



K, 



c 



0l±) 



6'(+) = 27r, al") = 



(4.29) 



Here, L± are the maximum and minimum values of the angular-momentum integration 
variable for the contribution of the circle orbits, Kq is their curvature (see Appendix [E|) . 

(a + l)(a-2) 



K, 



c 



2fLc 



(4.30) 



12 (v^ 

For a = 4; the period Tc [Eq. KW) ]. action Sc [Eq. K2^ ]. curvature Kc [Eq. ( lOOj) 
and stability factor Fmc [Eq. (12. 9p ] for the circle orbits are identical to those obtained in 



Ref. 



18| . We used also the properties of the Jacobians for transformations of the different 



coordinates, in particular, given by Eq. flESp . The time- reversibility of the Hamiltonian was 
similarly taken into account in Eq. (I4.24p . as explained above. Note that after applying the 
stationary phase conditions r* = r^ [Eq. 14.18] andp* = [Eq. (I4.2ip ]. the angular momentum 
L of the circular orbits as function of the r and Pr becomes the stationary point L* = Lc 
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at the boundary of the classically accessible phase space. Notice also that the asymptotic 
Maslov phase is defined traditionally in terms of the Maslov index (j\.l(j [Eq. fl4.13p ]. There 
is again the two components of the Maslov phase in the ISPM trace formula f l4.24p for the 
MC orbits. One of them is the asymptotic constant part f l4.13p independent of the energy. 
Another part is the argument of the complex amplitudes Aj^j^j [Eq. fl4.27p ]. that changes 
continuously through the bifurcation points. The total Maslov phase for the circle POs is 
given by the sum of these two contributions. It ensures a smooth transition of the trace 
formula fl4.24p for the contribution of the circle POs through the bifurcation points. 

In the asymptotic limit of the finite non-zero integration boundaries, L_ — t- — oo and 
0^ — )■ oo, i.e., far from any bifurcations Obif (also from the HO symmetry breaking at 
a = 2), the expression (14.271) tends to the amplitude of the Gutzwiller trace formula for 



13| . The simplest approximation is L+ = Lq, L^ = 0; and O^ = 0, 



isolated orbits [3|, 

Or = 2%, which correspond to the total physical phase space accessible for the classical 



motion. The factors va + 2 in front of the curvature Kc in Eq. ( 14.29^ appear because of 
the frequency ratio Eq. (12. 8p for uj^p/ur for the circle orbits for any parameter a > 2. Taking 
Eq. 04.271) asymptotically far from the bifurcations ctbif [Eq. (12. lip ], one has the amplitude 
of the Gutzwiller trace formula: 



l(0) IT7\ . _}_ T^ 

27r h ^Mc 



AZiE)^--^^. (4.31) 



In this limit, the asymptotic Maslov index crjy/^ and 0|^ in Eq. (14. 8 p are given by Eq. (I4.26p . 
For the opposite limit to the bifurcations {Fmc -^ 0, when a — )■ ctbif), one finds that 



the both arguments of the second error function in Eq. (I4.27P tend to zero as a/ | -Fmc I? see 
Eq. (14.291) . The Gutzwiller stability factor -Fa/c, going to zero, is exactly canceled by the 
same one in the denominator, and we arrive at 

Amc\E) -^- 



x^^i{K:iic^K%) ^"'''"- (4-32) 

Thus, in contrast to the Gutzwiller SSPM, one obtains the finite result at the bifurca- 
tions within the ISPM. Notice that the enhancement in order of fiT^I'^ with respect to the 
Gutzwiller asymptotic amplitude (I4.3ip takes place locally near the bifurcation points. Note 
also that at the circular billiard limit, when Kc — > C)0 (separatrix), one finds a continuous 
limit which is zero in the case of the PRL potential. 
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C. Total trace formula for the oscillating level density 



The total semiclassical oscillating (shell) correction to the level density (13. ip for the RPL 
potentials in two dimensions is thus given by 



,(1) 



(0), 



S9sci{E) = 6gl^{E) + 69i:>{E) 



(4.33) 



where 



«(i^) = ReJ]4'g(E) 



(/c). 



PO 



X exp 






-Spo(£/) -i-^cr^ 



dic) 



(4.34) 



i(^) 



The amplitudes Ap^ [see Eqs. (14. 9 p for /C = 1 and (I4.27P for /C = 0], actions S'po, Maslov 



indexes cTpQ, and constant phases 0)^ [Eqs. (I4.13P and (I4.26P 



were specified in the two 



previous sections IIVAI and IIVBI respectively. 

Using the scale invariance (12. 2p . one may factorize the action integral S-po{E), 



Spo{E) 



p ■ dr 



PO(E) 



E 
E'o 



2^Q 



p ■ dr 



PO{E=Eo) 



FT 



(4.35) 



In the last equation, we define the scaled energy e and scaled period t^q by 

p ■ dr . (4.36) 



E 
e'o 



-+- 

2^ a 



'po 



IPO{E=Eo) 

To realize the advantage of the scaling invariance (12.21) . it is helpful to use the scaled energy 
(period) in place of the corresponding original variables. For the HO one has a = 2, and 
the scaled energy and period are proportional to the unsealed quantities. For the cavity 
potential {a — )■ oo), they are proportional to the momentum p and length £po, respectively. 
Using the transformation of the energy E to the scaled energy e, one can introduce the 
dimensionless scaled-energy level density. The advantage of this transformation is that a 
nice plateau condition is always found in the Strutinsky SCM smoothing procedure by using 
the scaled spectrum Si (see Refs. ll|, |22| for the case of the billiard limit a — )■ oo). Then, 
one can use a simple relation between the original and scaled-energy level densities, 

dE 



g{e) = J2Ke-e.)=g{E) 



de 



(4.37) 



For the corresponding semiclassical formula for the oscillating part of the level density, one 
finds 



60, 



scl '' 



d^ 



^Sg^^\E) = j:6glSie) 



de 



Re Y^ ApQ (e) exp 

PO 



PO 



— FT 









l(/C) 



dE 



(ic). 



A^ciie) = -T-A'p^^iE). 



(4.38) 



^pov-. ^^ 

The simple form of the phase function fl4.35p enables us also to make easy use of the Fourier 
transformation technique. The Fourier transform of the semiclassical scaled-energy level 
density with respect to the scaled period r is given by 

F{t) = I de g{e)e''^"' ^ F,{t) + ^Ipo5(r - r^J , (4.39) 

•^ PO 

which exhibits peaks at periodic orbits r = TpQ. Fq{t) represents the Fourier transform of the 
smooth Thomas-Fermi (TF) level density and has a peak at r = related to the zero-action 
trajectory [1^. Thus, from the Fourier transform of the scaled-energy quantum-mechanical 
level density. 



Fir) 



E^(^ 



^ier/h d£ = ^ . 



,ieiT/h 



Eo 



2^ a 



(4.40) 



one can directly extract the information about classical PO contributions. The trace formula 
(I4.33P has the correct asymptotic SSPM limits to the Berry and Tabor results 04. 8p . (I4.14p 
for /C = 1 polygon-like (including the diameters) and to the Gutzwiller trace formula 04.241) . 
()4.3ip for /C = circle POs. As shown in the sections IIV Al and IIV Bt one obtains also the 



imit of the trace formula [Eqs. (I4.33P and 04.34p ] to that of the circular billiard a — t- oo 
13l . 1271 . In this limit one has obviously zero for the circular-orbit contributions as for the 



potential barrier separatrix in the IHH potential 



23| 



D. Harmonic oscillator limit 



In the isotropic harmonic oscillator limit [a — > 2 in the power-law potential 01.11) ]. the 
energy surface is simplified to the linear function in actions, 

E = U)rIr + U)vI^= UJ^ (2 Ir + L) . (4.41) 
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Therefore, in this hmit the curvature Kpo for all POs [including the maximum value L = 
Lc = E/uj^ for the circle orbits, Eq. f l4.30p . and L = for diameter ones, Eq. flD3l) ] and 
stability factor Fmc [Eq. ( 12. 9p ] for the MC orbits turn into zero. However, there is no 
singularities in the ISPM trace formulas fl4.8p for the contributions of all /C = 1 families and 
fl4.24p for the circle orbits in the limits K-po — )> and F^c ~^ 0. The arguments of both 
error functions, oc \/Kc and oc \J FmcI^c in Eq. fl4.27p . for instance, approach zero and 
singularities are canceled with the same ones in the denominators of the multipliers in front 
of them, and similarly, in Eq. (14. 8 P for one error function; see Eqs. (14. 9p . f l4.12p . f l4.27p and 
f l4.29p with the help of Eq. fl4.32p . Therefore, one has a continuous limit of the total trace 
formula fl4.33p for a — )■ 2. Moreover, in this limit, one obtains exactly the same half of the 
HO trace formula for the MC orbit contribution f l4.24p and the M(2, 1) diameter one [Eq. 
f l4.8p ] u p to the relatively small higher-order corrections in fi [see also Eq. (3.68) of Sec. 3.2.4 



in Ref . 



13|] 



.(0) ^p^ ^ ^A.(2)^p^ „(i) iT?\^^RrP) 



aTilicyiE) ^ 2 ^SnoiE), g\lj,^{E) ^ - 5g'^>o(E) . (4.42) 

Here, {MC} and {MD} represent sum of all repetitions of circle and diameter orbits, 
M = 1,2, ..., respectively. Thus, the HO limit of the sum of the circle and diameter orbit 
contributions is exactly analytically given by the HO trace formula. We point out that for 
a — > 2, the contributions of circle MC and diameter M(2, 1) orbits encounter local increases 
of the degeneracies /C by 2 and 1 units, respectively. 

As noted above, in the HO limit a — ;■ 2, only the diameter M(2, 1) and the circle MC 
(both with repetitions) survive, and they form /C = 2 families in the HO potential. Taking 
into account also that the angular momentum for the diameters is always zero, L* = 0, and 
for the circle orbits L* = Lq, we shall assume that the integration over L for the diameters 
is performed from L„ = to L+ = Lc/2 and for the circle orbits from L„ = to L+ = Lq, 



such that they give natura. 



ly equivalent contributions into the HO trace formula, as shown 



in Eq. (14. 42 p . see also Ref. 23|. The difference is in the integration limits for the circle orbits 
[see Eq. (I4.29p ]. in contrast to Eqs. (I4.12p and ( lD2p for the diameter boundaries. Notice that 
the contribution of the polygon-like one-parametric orbits, 5g^^\E), disappears in the ghost 
HO limit. Thus, one obtains the continuous transition of the oscillating part of the ISPM 
level density 6gsciiE) through all bifurcation points, including the HO symmetry breaking. 
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E. Averaged level density 

For comparison with the quantum level densities obtained by the SCM, we need also 
to average the trace formula fl4.33p over the spectrum. As this trace formula is given by 
the sum of the separating periodic-orbit terms everywhere (including the bifurcations), one 
can approximately take the integral of the local averaging over energies (folding integral) 
analytically in terms of the Gaussian weight factor at the Gaussian width parameter F <^ 
Ep. As result, for the averaged density 5g^{E) with this weight function, one obtains (see 



m 



Refs. J5|,ll2|,ll3|) 

5g,{E) = J2^9,oiE) exp [- (Ur/hf] . (4.43) 

PO 

The total level density is calculated analytically by adding the TF smooth component g^^ (E) 

m, 

g,{E)=g^^{E) + 5g,iE), (4.44) 



where 



grAE) = j^ J dv J dp 5 iE-Hip,r)) 



2 I f eV^" 

' ^ . (4.45) 



2!f 2Eo \Eo 

Here, H{p, r) = p^ /2m+V{r), p{r) is given by Eq. (12. 4p . and rmax is the turning point as one 
of solutions of the equation V{r) = E. For the RPL potential (12. ip . Tmax = Ro{E/ E^y/"", 
and we used the relation Eq = h'^/mRQ in the last equation of (14. 45 p . 

For the averaged scaled-energy level density, one obtains the semiclassical trace formula: 



1 



EE'^#o(^)^^p -( 



K.=0 PO 



'7po7^2 
2h 



(4.46) 



where 5Q^Q{e) is given by Eq. (I4.38p . 7 is the dimensionless Gaussian width parameter used 
for the coarse-graining over the scaled spectrum ei. For the Thomas- Fermi density, one 
obtains the simple expression: 

fe(.) = 3..mf = ^e. (4.47) 
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F. The shell correction energies 



The PO expansion for the semiclassical sheU correction energies SUsd is given by 



13 



El 



11 



(5f/sei = 2j]|^5ft,jE^), (4.48) 



PO ^PO 



where t^^ = MT-po{Ef) is the period of particle motion along the PO (taking into account 
its repetition number M) at the Fermi energy E = Ep. The Fermi energy Ep as function 
of the particle number A^ is determined by the particle number conservation, 

N = 2y^ni = 2 [ \e g{E) , (4.49) 

where rii = 6{Ef — Ei) are the occupation numbers. The factors 2 in Eqs. (14.481) and (I4.49p 
account for the spin degeneracy of Fermi particles with spin 1/2. 

Note that the shell correction energies 6U which are the observed physical quantities 
do not contain an arbitrary averaging parameter F, in contrast to the level density gp{E). 
The convergence of the PO sum (I4.48P is ensured by the additional factor in front of the 
oscillating density components 6gpQ which is inversely proportional to square of the PO 
period tpQ . This means that we need shorter orbits if they occupy enough large phase-space 
volume. 

In the quantum SCM calculations, the shell correction energies are usually obtained by 
extracting the oscillating part from a sum of the single-particle energies, 

U = 2^niEi. (4.50) 

i 

However, the direct application of the SCM average procedure to the spectra Ei of RPL 
potentials (except for the HO limit) does not give any good plateau condition as for the 
level density g{E) in Eq. (I4.37p . On the other hand, there is no problems in calculations of 
the shell-correction energies because one can always find rather a good plateau in the SCM 
application to a helpful sum of the single-particle scaled energies Ei, 

U = 2j2niei, (4.51) 

i 

which is simply related to the original shell correction energy. Indeed, applying exactly the 
same derivations of Eq. (I4.48P to the semiclassical trace formula for the oscillating part of 
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Eq. fl4.5ip . one gets 

5W,ci = 2j]4-^6;po(£j. (4.52) 

PO ^Po 

Here, the scaled Fermi energy e^ is determined by 

N = 2 [^ g{e)de . (4.53) 

Jo 

Using now the obvious relations 



dE/de' ^°' ' dE/de 

in Eq. (I4.48P one obtains 

SUsci = (^) SU,,i . (4.54) 

Thus, we arrive at the simple relation between the original shell- correct ion energy 6U [Eq. 
(I4.48P ] and the scaled one SU, valid for both semiclassical and quantum (neglecting the 
second order terms in the shell fluctuations of the Fermi energy) calculations: 

dE\ _, „ 2a (a-2)/{a+2). 



SU=i:^] 5U = Eo^^e'p^""^''5U. (4.55) 

This relation can be also directly obtained by using the standard quantum SCM relations 
of the flrst-order shell-correction energy 6U to the oscillating part of the level density Sg{E) 
up to the same second order terms in the Fermi energy oscillations [2], 

rEp 

SU = 2y] SuiEi = 2 dE {E- Ep) 5g{E) , (4.56) 

i Jo 

and corresponding ones for the scaled quantities, 

5W = 2 V 6n,ei = 2 f ^ de {e - e^) 5g{e) . (4.57) 

In these derivations, 5ni = Ui — rii, rii are the smooth occupation numbers of the SCM. We 
applied also the usual transformations from the Fermi energies to the particle numbers by 
using Eqs. (I4.49P and (I4.53p . as well as the deflnitions of the averaged Fermi energy Ep, and 
the scaled one ep, 

N = 2 f ^ dE g{E) = 2 f ^ de g{e) . (4.58) 

Jo Jo 

23 




6 8 10 12 14 16 rv 18 20 



a 



FIG. 2. Moduli of the scaled amplitudes |^po(e)| as functions of a for the primitive (M = 1) circle 
C [see Eq. (g^ZD], diameter (2, 1) and triangle-hke (3, 1) [see Eq. (gj])] POs, in units of e^/^ ^t the 
scaled energy e = 40. The panel (a) shows the HO limits (a — )■ 2) of the improved stationary-phase 
method (ISPM, solid), and standard stationary-phase method (SSPM, dashed curve) amplitudes 
for the circle C (thin) and diameter (2, 1) (thick curve) orbits. The filled circle denotes one half 
of the HO amplitude ([442]) at a = 2. In the panel (h), ISPM (sohd) and SSPM (dashed curve) 
amplitudes for the circle C (thin) and triangle-like (3, 1) (thick curve) POs are shown. 

V. AMPLITUDE ENHANCEMENT AND COMPARISON WITH QUANTUM 
RESULTS 



A remarkable enhancement of the ISPM amplitudes in PO sums for the oscillating level 
density (I4.34p and shell correction energy (14.481) is the general characteristic property of the 
periodic orbits, due to their bifurcation (symmetry-breaking) scenarios where the newborn 
POs emerge from the existing parent ones. In Fig. [2], the scaled amplitudes j^pol, divided 
by e^l'^ to normalize the energy dependence for /C = 1 orbits, are presented for several 
shortest POs as functions of the power parameter a in order to show the typical bifurcation 
enhancement phenomena. In Fig. [2](^ay), the enhancement of the primitive diameter (2, 1) 
amplitudes |^(2,i)| [Eq. (14. 9p ]. and those [.AmcI [Eq. (14.271) ] for the primitive circle orbit 
C are clearly seen in the HO limit a — )■ 2; see also Eq. (I4.38p . Figure W(h) shows the 
enhancement of the shortest orbit C around the bifurcation point a = 7, and the birth 
of the triangle-like orbit (3,1) there. Note that the ISPM amplitude |^(3,i)| [Eq. (14. 9p ] of 
the (3, 1) orbit keeps its magnitude up to rather a large value of a above the bifurcation 
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FIG. 3. The oscillating part of the level density 5Q^{£) in units of e^'"^ vs the scaled energy e for 
a = 4.0 (a), and a = 6.0 (b), at the (dimensionless) width parameter 7 = 0.1 in the Gaussian 
coarse-graining over the scaled energies: QM (solid curve) is the quantum-mechanical results (using 
the Strutinsky SCM with the averaging width 7 ~ 3 and curvature correction polynomial of the 
6th power); ISPM (dashed curve) is the semiclassical results using the improved stationary phase 
method. 

(a > ftbif)- The ISPM amplitude for the circle orbit C exhibits a remarkable enhancement 
at the bifurcation point a = 7. The divergence of its SSPM amplitude at the bifurcation 
point is successfully removed. As also seen from Fig. ^b), the ISPM amplitude for the 
(3, 1) PO is continuously changed through this bifurcation, in contrast to the discontinuity 
of the SSPM amplitude. This orbit exists, in fact, only at a > 7, and the amplitude in the 
region a < 7 is due to the formal stationary point which has no direct sense in the classical 
dynamics. Therefore, the corresponding PO is called usually as a ghost orbit 
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13 



23(. An 



oscillatory behavior of the amplitude |-4.(3^i)| in the ghost region far from the bifurcation has 
no physical significance, since it is washed out in the averaged coarse-graining level density 



by an rapidly oscillating phase of the complex amplitude ^(3,1); see Refs. [Ul, l23|. These 



ghost amplitude oscillations are suppressed even more b y u sing higher order expansions in 
the phase and amplitudes in a more precise ISPM (Ref. ll|). 

Figures [3H7] show the oscillating part of the semiclassical scaled-energy level density 6Q.y{e) 
[Eq. (I4.46P ] in units of e^^^ as functions of the scaled energy e for several values of the 
power parameter a and the averaging width 7. The ISPM semiclassical results show good 
agreement with the quantum mechanical (QM) ones for a transition from the gross to fine 
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FIG. 4. The same as in Fig. [3] for a = 4.0 but with other width parameters, 7 = 0.03 (a,h) and 
7 = 0.2 (c,d). Panels (a,c): QM (sohd hne) is the quantum-mechanical results; ISPM (dashed) is 
the semiclassical results using the improved stationary phase method; SSPM (dotted) is the results 
of the standard stationary phase method [shown only in the panel (a)]. Panels (h,d): {MD} 
(dashed) is the contribution of the diameters (including their repetitions) and {MP} (thin solid) 
for other /C = 1 polygon-like POs. 

resolutions of the spectra. The QM calculations are carried out by the use of the standard 
Strutinsky averaging over the scaled energy e, in which we find a good plateau around the 
Gaussian averaging width 7 = 2 — 3 with the even curvature correction polynomials of 4th 
to 8th powers. 

For the powers a = 4.0 and 6.0, one finds a good agreement with the SSPM asymptotic 
behavior [see Figs.|4]('aj and[5](^aj] because they are sufficiently far from the bifurcation points 
a = 4.25 and 7.0 which correspond to the birth of the star-like (5, 2) and triangle-like (3, 1) 
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FIG. 6. The same as in Fig. [5] for a = 4.25, but with the width parameters, 7 = 0.1 (a), and 
7 = 0.2 (b). 



POs (see Figs. [6] and [7] below for these bifurcations). For the gross shell structure (7 ~ 0.2 
at a = 4.0), only the shortest orbits (mainly a few shortest diameters) give the leading 
contributions. (This is in contrast to the 3D case where the circular orbits become also 
important; see Refs. |l3l. llSl].) For instance, the gross shell structure in terms of the shortest 
POs for a = 6.0 — 7.0 manifests at larger 7 > 0.3, unlike for the powers a = 4.0 — 4.25. With 
decreasing the averaging parameter 7 and increasing the power value a, the POs for larger 
scaled periods r [or actions S, see Eq. (I4.35P ] become more significant [compare Figs. ^b,d) 
SindU](b,d)]. In the case of the fine shell structure (e.g., 7 ~ 0.03) the dominant contributions 
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FIG. 7. The same as in Fig. U but for a = 7.0 and other width parameters, 7 = 0.1 (a,b) and 
7 = 0.2 (c,d). 

are due to the bifurcating 2D /C = 1 POs [polygon-like POs denoted byj MP}; see Fig.|4](^&^]. 
(This is similar to the situations in the elliptic 22| and spheroidal ll| cavities, and in the 
IHH potential |23|.) However, the interference of much longer these one-parametric POs 
[such as M(7, 3) for a = 4.0 or M(5, 2) for a = 6.0] with a lot of the M(2, 1) diameters 
explain some peaks, too. For smaller a = 4.0 and 4.25, the circle orbit contributions are 
not shown because they are not significant at these diffuseness parameters in the 2D case. 
(This is different from the 3D case; see Ref. 



for the trace formulas based on the uniform 



approximation using the classical perturbation approach.) These contributions into the trace 
formula (I4.46P are increasing functions of a, and they become significant at a > 7 even for 
the 2D case [see Fig. [7](^&,(ij]. An intermediate situation between the gross- and fine- shell 
structures, where all of POs become significant, are shown too at 7 = 0.1 in Figs. [3]and[6l 
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and in Fig. U^b,d) at 7 = 0.1 and 0.2. Note that our full analytical expressions (accessible 
for any long periodic orbits) for the classical PO characteristics at a = 4 and 6 are quite 
useful in the simple ISPM calculations of the oscillating level density with a good accuracy 
up to the fine spectrum-structure resolutions by using, for instance, 7 ~ 0.03 and 0.1. 

Figures Eland [7] show a nice agreement of the fine-resolved semiclassical and quantum level 
densities 6Q^{e) as functions of the scaled energy e at the critical bifurcation points a = 4.25 
and 7.0 for the births of the star-like (5, 2) and triangle-like (3, 1) orbits, respectively. 

Figures [8] and [9] show the scaled shell correction energies 6U [Eqs. fl4.52p for the semi- 
classical and fl4.57p for the quantum results], normalized by the factor Sp , as functions 
of the particle number variable A^^/^. A good plateau is realized for the QM calculations 
of the scaled shell-correction energies [see the first equation in Eq. fl4.57p ] near the same 
averaging parameters 7 and curvature corrections as mentioned above. In the semiclassical 
calculations, the Fermi level e^ is determined by the particle number conservation fl4.53p 
with using the coarse-grained scaled-energy POT level density. 
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The oscillating ISPM components 6Q^ idi^) were determined above by Eqs. (14.461) and (14.381) . 
We evaluated the Fermi level e^ (N) by varying the averaging width 7 and found that there 
is no essential sensitivity within the interval of smaller 7 (7 ;^ 0.1 — 0.2). Moreover, even 
the TF density ^tf(£^) [Eq. (I4.58J in Eq. (I4.53P with Q{e) ~ ^tf(£^) provides us a good 
value of Ep in the POT calculations of the shell correction energies (I4.52p . The PO sums 
at a = 7.0 converge for the shell correction density (14.461) by using the averaging width 
7 = 0.2 of a fine shell-structure resolution, and for the shell correction energies (14.521) with 
taking into account the same major simplest POs [about 4 repetition numbers (M = 4) 
for the circle and diameter orbits, and a few first simplest other /C = 1 (P) POs, such as 
(3,1), (5,2), (7,3) and (8,3); compare Figs. [9]^c,d) with [71^c,d)]. For smaller diffuseness, 
4 < a < 6, one has a similar PO convergence relation with the same 7 ~ 0.2, but with much 
smaller contributions of the circular orbits. However, the dominating (/C = 1) PO families 
(P) are the (5, 2), (7, 3) and (7,3) POs at a = 4.25 — 6.0 and 4.0, respectively [see, for 
example. Figs. [H!^a,b) and|l](c,d)]. As seen from Figs. [6], Hand [9], we obtain a nice agreement 
between the semiclassical (ISPM, dashed) and quantum (QM, solid curve) results exactly 
at the bifurcations a = 4.25 and 7.0. Notice that the dominating contributions in these 
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FIG. 8. Scaled shell correction energies 6V(, normalized by the factor £p , as functions of square 
root of the particle number N^'"^ at the values of a, where the full analytical formulas are obtained 
for a = 4.0 (a,b) and a = 6.0 (c,d). Panels (a,c): QM (solid curve) represents the quantum- 
mechanical results using the Strutinsky SCM, and ISPM (dashed curve) shows the semiclassical 
result using the TF approximation in the calculation of N{eF) by Eqs. (|4.53|) and (j5.ip . Panels 
(h,d): the contributions of several POs into the shell correction energy 5U are shown. Other 
notations are the same as in Figs. [Hand [71 

semiclassical results at the bifurcation point a = 7.0 are coming from the interference of 
the bifurcating circle C and newborn (3, 1) orbits with the simplest diameters. As shown 
in Figs. U^d) and W^d), one can see more pronouncedly that the circle C and triangle-like 
(3, 1) orbits are mainly in phase, but the diameter (2, 1) is sometimes in phase to them 
and sometimes out of phase. Thus, the occurrence of a characteristic beating pattern in 
the level density amplitude at a = 7.0 is due to the interference of the bifurcating orbits 
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FIG. 9. The same as in Fig. [HJ a = 4.25 for the (5,2) bifurcation, and a = 7.0 for the (3,1) 
bifurcation are shown in panels (a,b) and (c,d), respectively. 

C and (3, 1) with the shortest diameter (2, 1) having all the amplitude of the same order 
in magnitude but different phases. The bifurcating circle 2C and star-like (5, 2) orbits [as 
expected from the enhancement of the amplitudes of the circular C and triangular-like (3, 1) 
POs in Fig. [2] are more important for a = 4.25, though the primitive diameters become 
significant much compared to the bifurcation case a = 7.0. The POs (3, 1) and (5,2) yield 
more contributions near their bifurcation values of a, and even more on the right-hand side 
(a > ftbif) in a wide region of a as mentioned above. The partner bifurcation orbits {C, 
(3, 1)} and {2C, (5, 2)}, taken together with the simple diameter (2, 1), give essential ISPM 
contributions of about the same order of magnitude and phase in Figs. El [3 and IHl as seen 
for example in Figs. ^b,d) and^d) for the same a = 7.0. The diameter ISPM contributions 
are close to the SSPM asymptotic ones near the bifurcation points a = 7.0 and 4.25 (as for 
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FIG. 10. Moduli of the Fourier transform |-F(t)| of the quantum scaled-energy level density (|4.40p 
as functions of the dimensionless variable r are plotted for several values of a; MC and M{nr, n^) 
indicate the classical POs corresponding to each peak (see Fig. [1]). 



a = 4.0 and 6.0) because they are sufficiently far from their single symmetry-breaking point 
at the harmonic oscillator value a = 2. 

Figure [lO] shows the Fourier transform of the quantum-mechanical scaled-energy level 
density, see Eq. f l4.40p . For a smaller a = 2.1, the diameter (2, 1) orbit gives the dominant 
contribution to the gross-shell structure as the shortest POs; see the peak at r ~ 5.0. With 
increasing a, the amplitude of the circle orbit becomes again larger due to a prominent 
enhancement around the bifurcation point (r ~ 6.2 at ctbif = 7.0). Notice that the newborn 
POs (3, 1), (5,2), (7,3) and (8,3) give comparable contributions at a = 7.0 [similarly, (5,2) 
and (7, 3) for the bifurcation a = 4.25] in nice agreement with the quantum Fourier spectra 
in Fig. [TOl The contributions of the newborn triangle-like orbit family (3, 1) having relatively 
a smaller scaled period t and higher degeneracy /C = 1 become important and dominating 
for larger a > abif = 7. The newborn (3, 1) peak cannot be distinguished from the parent 
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circle C orbit near the bifurcation point abif as well as the diameter and circle orbits at a 



close to the HO limit, see Ref. 12|. We emphasize that the shell correction energies 5U are 
similar to the oscillating parts of the level density averaged over the spectrum by using the 
Gaussian width 7 = 0.2 at a = 4.0 and 4.25, which have mainly the gross-shell structure 
due to the shortest diameters. However, for this 7 the fine-resolved shell structures (due to 
their interference with the other polygon-like and circular POs) are pronounced at larger 
powers near a = 6.0, and especially, 7.0. 

Notice also that we do not show the numerical comparison of the ISPM trace formulas 
fl4.46p and fl4.52p vs quantum results for the HO limit a — )■ 2 because, as shown analytically 
around E g. (I4.42p . these formulas exactly analytically approach the corresponding HO trace 



formulas 13|. 



VI. CONCLUSIONS 

We presented a semiclassical theory of quantum oscillations of the level density and energy 
shell corrections for a class of radial power-law potentials which turn out as good approx- 
imations to the realistic Woods-Saxon potential in the spatial region where the particles 
are bound. The advantage of the RPL potentials is that, in spite of its diffuse surface, the 
classical dynamics scale with simple powers of the energy, which makes it particularly easy 
for POT calculations, sometimes even completely analytically for a = 4 and 6, beyond the 
well known HO case a = 2. The quantum Fourier spectra yield directly the contributions 
of the leading classical POs with the specific periods and actions into the trace formulas. 

We described the main PO properties of the classical dynamics in the RPL potentials as 
the key quantities of the POT. Taking the simplest two-dimensional RPL Hamiltonian we 
developed the semiclassical trace formulae for any its power a, and studied various limits of 
a (the harmonic oscillator potential for a = 2 and the cavity potential for a — > 00). The 
completely analytical results were obtained for the RPL powers a = 4 and 6. This can be 
applied for both 2D and 3D cases and allow us to far-going fine-resolved shell structures at 
7 = 0.03 — 0.1. This POT is based upon extended Gutzwiller's trace formula, that connects 
the level density of a quantum system to a sum over POs of the corresponding classical 
system. It was applied to express the shell correction energy 6U of a finite fermion system 
in terms of POs. We obtained good agreement between the ISPM semiclassical and quantum- 
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mechanical results for the level densities and energy shell corrections at several critical powers 
of the RPL potentials. For the powers a = 4 and 6, we found also good agreement of the 
ISPM trace formulas with the SSPM ones. The strong amplitude-enhancement phenomena 
at the bifurcation points a = 7 and 4.25 in the oscillating (shell) components of the level 
density and energy were observed in the remarkable agreement with the peaks of the Fourier 
spectra. We found a significant influence of the PO bifurcations on the main characteristics 
(oscillating components of the level densities and energy-shell corrections) of a fermionic 
quantum system. They leave signatures in its energy spectrum (visualized, e.g., by its 
Fourier transform), and hence, its shell structure. We have presented a general method 
to incorporate bifurcations in the POT, employing the ISPM based on the catastrophe 
theory of Fedoryuk and Maslov, and hereby, overcoming the divergence of the semiclassical 
amplitudes of the Gutzwiller theory and their discontinuity in the Berry& Tabor approach at 
bifurcations. The improved semiclassical amplitudes typically exhibit a clear enhancement 
near a bifurcation and on right side of it, where new orbits emerge, which is of the order ^~^'^ 
in the semiclassical parameter h. This, in turn, leads to the enhanced shell structure effects. 
Bifurcations are treated, again, in the ISPM leading to the semiclassical enhancement of 
the orbit amplitudes. The trace formulae are shown numerically to give good agreement 
with the quantum-mechanical level density oscillations for the gross- (larger coarse-graining 
averaging parameter 7 and a few shortest POs), and the fine- resolved (smaller 7 and longer 
bifurcating POs) shell structures, also for the shell- correct ion energies independent of 7. We 
found similar gross shell structures in the shell-correction energies with the corresponding 
averaged densities at 7 ~ 0.2 for smaller powers a = 4 - 4.25 of the gross shell structure, and 
for larger a > 6 of the fine shell structure. The fine-resolved and gross shell structures were 
found in the oscillating averaged densities at smaller averaging parameters 7 = 0.03 - 0.1 
and at larger 7 > 0.2 - 0.3, respectively. The fine-resolved shell structure for larger powers, 
a > 6, occurs in a larger interval, 7 = 0.003 - 0.2, including the essential contributions 
of the circle orbits along with the polygon-like and diameter orbits. Full explicit analytical 
expressions for the diameters and circle orbit contributions into the trace formula as functions 
of the diffuseness potential parameter a are specified too. 

We may state that the semiclassical POT is well capable of explaining the main features 
of quantum shell structure in terms of a few classical periodic orbits for the RPL potentials. 
Bifurcations of POs are not simply an obstacle of the semiclassical theory, but they leave 
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clear signatures both in the quantum Fourier spectra and the locations of minima of the 
shell-correction energy 6U which are plotted versus particle number and powers. 

For prospectives, we intend a further study of shell structures in the 3D RPL potentials, 
applying both the ISPM and uniform approximations to treat the bifurcations, by vary- 
ing continuously the power parameter a from 2 (harmonic oscillator) to oo (3D spherical 
billiards). 
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Appendix A: The stability factor, bifurcation powers and frequencies 

Let us consider in more details the non-linear classical dynamics in the RPL Hamiltonian 
(12. ip for any real a > 2. The critical values of the radial coordinate r = r^ and angular 
momentum L = Lq for the circle orbit (C) are determined by the solutions of the system 
of the two equations with respect to r and L : 

J^{r, L) = 0, ^^ = 0, where J^(r, L) = pl{r, L) , (Al) 

or 

see Eq. (12. 4p . In the internal region where the stable orbits in the radial direction exist, 
one has a nonzero J-"^ = d'^J^ {r^,Lc) /dr"^ = —AmEa/r^ < (VJ > 0). First equation in 
Eq. ( JAip means that there is no radial velocity, f = 0, and the next equation is that the 
radial force is equilibrating by the centrifugal force. For instance, for the potential fll.lD . 



the solutions of the two these equations are the radius r^ and angular momentum Lc 15| . 

f 2E \ ^/" 
'■'^ = ^° i(2TT^j • i- = Pfc)r„. (A2) 

Using Eqs. (12.70 at L = L^ and flA2p for r^ and Lc, one finds [l 



c.,=c.,(L = Lc) = ,/^(^^^^) . (A3) 

Applying now the second order expansion in r — r^ to Eq. (12. 4p . one gets the first-order 
ordinary differential equation, which determines the radial CT r{t) locally near the circle 
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PO r = r^: 



r = ± 



-nil 



C ) 1 



7'^ = -2m 



r\f-^ \ m ■r'-^ 



2mr' 



(A4) 



r=r^,L=Lc 

Integrating the dynamical equation in Eq. flA4p . one obtains 



r(t) = Tp + (r' — r^^) exp ± 



T" 

2 777,2 



t 



y = r (t = t' = 0)] . (A5) 

In the stable case, J-"^ < in Eq. (lASp for the CT r(t) locally near the circle orbit r = r^^ 
one writes 

r(t) = r^ + (r' - r^) exp {±iQc t) , (A6) 

where i?^ is a positive radial frequency Ur a,t L = Lc [see Eq. fl2.7p ]. 



a 



c 



T^/{2m^ 



Ur(L 



^C) 



(A7) 



For the potential ( II. ip . this quantity is given by |l5 | 

f2c = 



'a(2 + a)Eo 



mRl 



2E 



>0 



[{2 + a)Eo\ 

From Eq. (]A6P after the period T^ along the primitive circle orbit, 

2TTn,„ 



T, 



c 



t" -t' = t" 



UJ^ 



■, ^c 



u^{L = Lc) 



one finds 

6r" = r" ~ r^ = Sr' exp {±if2c Tc) , Sr' = r' — r^ . 

The eigenvalues of the stability matrix Aic for M = 1 in Eq. (12.90 are given by 

'dp':: 



dr" 
dr' 



exp{if2cTc) 



dp'j. 



exp {—iQc Ti 



c) 



(A8) 



(A9) 



(AlO) 



(All) 



p;, \ - i r / r' 

These two eigenvalues of the stability matrix are complex conjugated, as it must be, ac- 
cording to its general properties. So, as Qc is real [Qc > 0, according to Eqs. (1A7P and 
(lASp ] the circle orbit is isolated stable PO. Substituting the expressions (lAlip into the first 
equation in Eq. (12. 9p and using Eqs. ( 1A9P for the period Tc, (1A3P and (lASp for the C-orbit 
frequencies u^ and i7c, relatively, one obtains the last equation in Eq. (12. 9 p for the stability 
factor Fmc- 
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Appendix B: Scaling properties 

For convenience, let us consider the classical dynamics in terms of the variables in di- 
mensionless units m = Rq = Eq = 1. Due to the scaling property fl2.2p for the classical 
dynamics in the Hamiltonian f l2.ip . the energy dependence of the action /^ [Eq. fl2.5p ]. the 
angular momentum L, the frequency Ur [Eq. (12. 7p ] and the curvature K [Eq. f l2.12p ] can be 
expressed in terms of the simple powers of the scaled energy e, 

In particular, one can express these classical quantities through their values at e = 1 [E = 1), 

/, = /,(!)£, L = L{l)e, a;-^ = a;7^(l)£(2-")/(2+"), 

K = K{l)/e, (B2) 

where /«(!), -^(1), <-^r^{^), and K{1) are the corresponding quantities of Eq. ( IB2p . which are 
taken at "the scaled energy" e = 1. Therefore, due to the scaling properties (12. 2p and (IB2p . 
we need to calculate these classical dynamical quantities only at one value of the energy 
e = 1. For simplicity of notations, we shall omit the argument e = 1 everywhere, if it is not 
lead to misunderstandings. 

The energy surface Ir{L,E) for the potential ( 11.11) can be expressed explicitly in terms 
of the frequencies Uj. and u^^ [see left identity in Eq. (12. 8p ]. 

Ir = ^co;'-LfiL). (B3) 

a + 2 

To prove this identity, we express Eq. (12. 7p for u^^ in terms of the determinant, 

_i ^ fdiA ^ d{Ir,L) ^dh_dhdL_ 
"^^ \dEj^ d{E,L) dE OLdE' ^ ' 



Calculating directly the derivatives in this equation by using Eq. f lB2p . one obtains the 
expression for uj~^{V). Solving then this equation with respect to /r(l), one arrives at 
Eq. ( 1B3p . Differentiating the identity (1B3P term by term over L and using the definition for 
the ratio of frequencies f{L) in Eq. (12. Sp . for the curvature (I2.12p one finally obtains 

2a doj-^ 



K 



{a + 2)L dL 

a dTr _ 27r 



7r(a + 2)L dL ' w 
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T, = — , (B5) 



r 



where u^ ^ is obviously simpler quantity to differentiate over L than f{L), according to Eqs. 
(K7h and ([21]) with the help of the textbook Ref. ^, 

^;' = -^ / , , (B6) 

27rV2 ^„,i„ y^Q{x,L,a) 

g(x,L,a) = (l-x"/2) X-LV2, (B7) 

and X = r^. The turning points x^in{L,a) and Xma.xiL,a) are determined by the following 
algebraic equation: 

Q{x,L,a) = 0. (B8) 

Thus, we may calculate u~^ and f{L), and then, use Eqs. (1B3P and (IBSp for the radial action 
Ir, and curvature K at the scaled energy e = 1. Then, one obtains their energy dependence 
through the scaling equations flBl|) and flB2p . respectively. 

Appendix C: Full analytical classical dynamics for po^vers 4 and 6 



For the powers a = 4 and 6, the roots of function (]B7|) . in particular, the turning points 
Xmin and Xmax cau be obtained explicitly analytically. Therefore, one can find the explicit 
analytical expressions for the key quantities of the classical dynamics for the POT, namely, 
the radial frequency Ur [Eq. (12.71) ] (or the radial period T^), and the frequency ratio f{L) 



[Eq. (12.81) ] in terms of the elliptic integrals from Ref. 29|] (all in dimensionless units). Then, 
we calculate analytically also the curvature K through Eq. (1B5|) and action by Eq. (]B3|) . For 
a = 4, one has the cubic polynomial equation Q{x, L,a) = x — x^ — L^/2 = [Eqs. (IB7p and 
(jBSp ] for the three roots Xmm, a^max and xi which are given by the Cardano formulas explicitly 
as functions of L in the physical region L < Lq, ri < < rmin < rmax; Xq = r^. For the 
radial period T^ [Eqs. (IB 5 1) and (12. 7p ]. one obtains the analytical expression through these 



roots in terms of the complete elliptic integral F(-7r/2, k) of the first kind (see Refs. 18 



2a|) 



r, = H^ = ^^F(^,K). (ci) 



where 



V v-^max •^minj / l^a/jnax Xi) . y^ZJ 

For the ratio frequencies f{L) [Eq. (12. 8p ]. one finds 

/(L) = -y= ^ n f ""'-""- , ^ , (C3) 

TTV 2 XmaxV-^max — Xi \ ^max / 
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where n(n, k) is the complete elhptic integral of the 3rd kind 29 1. 



For a = 6, one has the polynomial equation of the 4th power, Q{x, L, 6) 



x — x 



-L^/2 



0, having the 4 roots [two complex conjugated Xi + 1x2 and xi — ix2, and again, two real 
positive roots, Xmm and x^aax', see Eqs. flB7p and fIBSp ]. The radial period T^ is determined 



through these roots by the expression [similar to Eq. (IClj . see Refs. 18|, |29 | 



where 



T 

J. r 



V2 



An (Xinax ~ Xi^ 



F.f.. 



\ ^ m a Y ^ Ti 



{A - Bf 



AAB 



(C4) 



(C5) 



A 



Xr, 



Xi 



X2, D — Y [Xxmn Xl) ~r X'. 



(C6) 



(We reduced the 4-power polynomial equation to a cubic one and obtained its 4 analytically 
given roots, mentioned above, in the explicit Cardano's form as functions of L). For f{L) 
of Eq. (12. 8 p at = 6, one obtains [29 1 



fiL) 



V2L (A + B) 



71 y AB (AXmin — ijXinax 



^ V2 / 2 1-/32 I ' 1 



/3^ 



/?■ 



:,«; 



(C7) 



where n((y9,?7,, k) is incomplete elliptic integral of the 3rd kind. 



/3 






/3i 



A-B 

A + B 



(C8) 



The curvatures K [Eq. (]B5] for a = 4 and 6 are determined by taking analytically the 
derivative of the radial period Tj. [Eqs. fICip and flC4p ] over L through the derivatives of the 
roots Xmin{L), x^s.^{L), Xi{L) and X2{L) by using Eqs. (lC2p . (ICSp . ( IC6D and f ICSD . and the 



explicit expression for the derivative of F(7r/2, n) over k (see Ref. 29||). The expressions for 



the curvatures K at the both powers a = 4 and 6 can be found in the closed analytical form 
through a rather bulky formulas, which contain the complete elliptic integrals of the 1st and 
2nd kind (for instance, by using Mathematica) . 
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Appendix D: Classical dynamics and boundaries for the diameters 

For the primitive diameter D = (2, 1), the action Sd (all in this appendix in dimensionless 
units) is specified analytically through the scaled period r^ and energy e by 



S 



D 



TdS , 



a + 2 \a 



^ ,1 1 

r(- + - 

2 a 



(Dl) 



where T{x) is the F-function of a real positive argument x. For the diameter PO boundaries, 
one can use the same L_ = 0, but L^ = bj^Lc, where 



"^o = 1 - - exp 



rHO T 



'D 



lA 



D 



A 



D 



^J-nMurKr) 



(D2) 



(see Ref. 23[ and more details in relation to the HO limit in Sec. IIVDI) . Lo = is the 
stationary point, L^^ = Lc/2 = £:/(2v2) is the upper angular momentum L+ for the D 
orbits in the limit a — )■ 2, in which 6^ — t- 1/2. In the semiclassical limit e ^ 1, one has 
6/5 —7- 1. Au is the Gaussian width of the transition region between these two asymptotic 
limits. The diameter-orbit curvature Kd for a > 2 at the stationary point L = L^ is given 
by 



K 



D 



F(l-l/a) 



eV2TTT{l/2-l/a) 



(D3) 



This exact analytical expression for the curvature K^, at any a was derived by using a power 
expansion in equation (JBSP for the turning points, and in the identity (1B3P over the variable 
proportional to L^ near L = up to the terms linear in L^. The Maslov phase for the 
diameter orbit was determined by Eq. (14.131) at n^ = 2 and n^ = 1. Note that for the limit 
a — )■ 4, the general expressions for the period r^ and action So [Eq. flDip ]. the Maslov index 
(Td [Eq. (14.131) ] with the same asymptotic (SSPM) limit of the constant part of the phase 
0|^ = —n{Ab£) — 3)/4 — )> — 7r/4, and the curvature K^) [Eq. (1D3P ] for the diameter (2, 1) are 



identical to those obtained earlier in Ref. 
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Appendix E: The boundaries and curvature for circle orbits 

For the arguments Zj, and Zi- of the error functions in Eq. fl4.27p . one originally has 






4Uc=V-i^Mc(-^^^--o), (El) 

where p* = is the stationary point, pr and r^^-* are maximal and minimal classically 
accessible values of pr and r as the finite integration limits for the corresponding variables. 
To express the integration boundaries (lEip in an invariant form through the curvature 
Kc ( 14.30p . and stability factor Fmc (ES]), one may use now the simple standard Jacobian 
transformations, and the definition of the angle variable 0'^ as canonically conjugated one 
with respect to the radial action variable Ir by means of the corresponding generating 
function. In these transformations, we apply simple linear relations: 

f dr'Y 
-'--*=(^j i^'r-O:). (E2) 

where we immediately recognize the Jacobian coefficients. Note that there is no cross- 
ing terms due to the isolated stationary point /* = 0, 6** = and to equations for the 
canonical transformations. At the stationary point for the isolated circle PO, one has 



f[Lc) = -{dIr/dL)L=Lc = -l/Vtt + 2; see Eqs. I^M) and ( KWt . For the transforma- 
tion of the derivative dr" /dO'^., one can apply the Liouville conservation of the phase space 
volume for the canonical variables to arrive at dr" /dO'^ = {dir / dL) / {dp'j. / dL) . We use also 
the following identities: 

Fmc = -JMC-^Mcl-^MciPr-^Vr) ) 

dPr ^ dp'^/dL ^ r%_ ^ ^ 
dp[. dpr'jdL r"p'l 

i.e., |i7cT(Pr7Pr)l = 1 at the PO conditions r' — )■ r" — )■ r^, pj. — ;■ p'^ — )■ 0. After such simple 
algebraic transformations, one obtains Eq. f l4.29p for the arguments of the error functions in 
Eq. (1127D- 

The expression (14.301) for the circle-orbit curvature Kc (in dimensionless units at e = 1) 
was obtained from expansion of f{L) [see Eq. (12. 8p ] as function of L in powers of Lq — L = e^ 
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up to the 2nd order terms in e. For this purpose, we have to solve first Eq. (IBSP for the 
turning points rmm and Tmax, which are the integration hmits in Eq. (12. 8p by using standard 
perturbation theory looking for the solutions in the following general form (r is taken below 
in units of Rq), 



re + cie + c^e + cse + c^e + 



r-min = rc- cie + C2^ - c-^e" + C4e^ + 



(E4) 



Existence of such form of the solutions follows from a symmetry of the equation flB8l) with 
respect to the change of the sign of e. Substituting these solutions into Eq. (IBSP for arbitrarily 
small e, one gets the system of the recurrent equations for the coefficients c^. The solutions 
of this system up to the 4th order in a perturbation parameter e is given by 



C\ 



Lc 

a 



C2 



a + 1 



■ Ci , C3 
„2 



6 ^^' ^-^ 72 



a-2)(2a + 5) 3 



C4 



(a + l)(4a2 + 8a + 13) 
1080 



'1 > 



(E5) 



and so on. We transform now the integration variable r in the integral of Eq. (12. 8 P for f{L) 
to y, r = r^(l — y), so that 



fiL) 



L 



n 



ymax 



dy 



y) ^/Q{y,L,a) 



where Q{y,L,a) is given by Eq. (IB 71) . 



(E6) 



Q{y,L,a) 



2rl 



1- 



a + 2 

iVmi^K - y)iy - yminWy) 



l-yY-Ll + 2Lce' - e' 



(E7) 



Here, 



?/max = Cie - C2e^ + 036^ - 046^, 

l/min = -Cie - Cae^ - c^e^ - c^e^ 



(E8) 



where Cn = Cn/r^. We use the last representation in Eq. f lE7p . introducing a new function 
'R-{y) of the new variable y to separate the singularities of the integrand in Eq. flE6p due to 
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the turning points. This integrand has to be integrated exactly by using a smooth function 
TZ{y) of y, which can be expanded in y at y = up to the second order, 

niy) = 7^(0) + 7^'(0)y + ln"iO)y' + ■■■. (E9) 

In order to get analytically the final result, we note that y in this expansion is of the order of 
e, according to Eq. flE8|) . Substituting then these expansions flE8l) and flE9P into very right 
of Eq. ( 1E7I) . we expand their middle in y at y = up to the 4th order. After the cancellation 
of e^ from both sides, and simple algebraic transformations, one has 

7^(0) = ^[l + e'c2{l- h)] + 0{e'), h = 3, 

'^ = 2h + 0{e'), ^ = 2h{3h + l) + 0{e'). (ElO) 

For the calculation of the circle orbit curvature Kq, we obviously need only quadratic terms 
in e [linear in {Lc — L)]. Therefore, one may neglect the e^ corrections in the second 
and third lines of Eq. flElOp because they are multiplied by y ~ e and y"^ ~ e^ in the 
expansion f lE9p . respectively. Substituting now expansions f lESP and f lEQP into the integral 
over y in Eq. ( 1E6|) . and taking 7^.(0) off the integral, one then expands to the second order 
all quantities of the integrand in y ~ e, except for (t/max — y){y — ymm) under the square 
root (in the denominator) which can be integrated exactly. Taking remaining integrals as 



/ o??/?/''/A/(ymax -y){y- ymhi) from yrnin to y^ax [see Eq. ([EH])], and then, expanding finally 
f{L) [Eq. (1E6P ] in e, we find that the linear terms exactly disappear. It must be the case 
because f{L) is an even function of e. Thus, the coefficient in front of e^ with the expressions 
for c„ (n = 1, 2, 3) from Eq. ( 1E5I) is just Eq. fl4.30p for the curvature Kc- We can also use this 
perturbation method for the calculation of the next order curvatures, for instance, the third 
derivative of Ir over L, which appears in expansion of the phase integral in the exponent up 
io the third order terms near the stationary points within a more precise (3rd-order) ISPM 
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